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INTRODUCTION 

A new method for estimating the power spectrum Ci from cosmic microwave 
background (CMB) maps was recently presented by the author^ and applied to the 
2 year COBE data, giving the results in Eigure 1. It was found that the spectral 
resolution for COBE could be more than doubled sX I = 15, thereby revealing 
previously unresolved features in the power spectrum. Whereas that paper was 
rather technical, the present paper shows that all qualitative features of this method 
can be understood from a simple analogy with quantum mechanics. 

There has been a surge of interest in the cosmic microwave background radiation 
(CMB) since the first anisotropies of assumed cosmological origin were detected by 
the COBE DMR experiment^. On the experimental front, scores of new experiments 
have been carried out and many more are planned or proposed for the near future. 

On the theoretical front, considerable progress has been made in understanding 
how the CMB power spectrum Ci depends on various cosmological model parame¬ 
ters, both analytically^ and quantitatively^. It is therefore quite timely to further 
strengthen the link between these two fronts, by better understanding how to ex¬ 
tract more accurate power spectra from experimental data. Early work on this 
problem®’® has recently been extended and applied to the 2 year COBE data^. 
When estimating power spectra, it is customary to place both vertical and hor¬ 
izontal error bars on the data points, as in Eigure 1. The former represent the 
uncertainty due to noise and cosmic variance, and the latter reflect the fact that 

Figure 1. The power spectrum 
observed by COBE. 

The observed multipoles are plot¬ 
ted with 1 — (7 vertical error bars, 
including both pixel noise and cos¬ 
mic variance. The horizontal er¬ 
ror bars show the r.m.s. window 
function widths. If the true power 
spectrum is given by n = 1 and 
Qrms,ps = 20yuK (heavy horizontal 
line), then the shaded region gives 
the 1 — (7 error bars and the double- 
shaded region shows the contribu¬ 
tion from cosmic variance. 

0 5 10 15 20 25 30 

1 



1 















an estimate of Ci inadvertently also receives contributions from other multipole 
moments. As is well-known, this unavoidable elfect is caused by incomplete sky 
coverage, which destroys the orthogonality of the spherical harmonics. For typical 
ground- and balloon-based experiments probing degree scales, the spectral blurring 
/S.l/1 tends to be of order unity, which makes it difficult to resolve details such 
as the number of Doppler peaks. By expanding the COBE sky map in spherical 
harmonics^, the spectral resolution /S.l/1 can be brought down to the order of 25%. 
The optimal method^ reduces the horizontal error bars still further, down to their 
theoretical minimum, which is seen to be ~ 1 for the COBE data in Eigure 1. 
It involves solving a generalized eigenvalue problem numerically, were the matrices 
involved depend on the survey geometry and the pixel variances. 


THE QUANTUM ANALOGY 


Although all the elements of realism included in the numerical technique^ are 
important when applying it to real data, one can in fact understand all of the 
qualitative features of what is going on by ignoring most of these complications. 

Eirst of all, let us ignore the fact that a real CMB sky map is pixelized. We let the 
function a;(r) denote the temperature fluctuation AT/T in the direction of the unit 
vector r. In the discrete real-world case, each pixel has some r.m.s. noise variance 
(7^ and effectively covers some solid angle Q of the sky. Eor our continuous analogue, 
let us define the noise variance function U(r) by U = cr^Q, where a and Q refer to a 
pixel in the direction r. If we neglect contamination problems, it is easy to see that 
V (?) is simply proportional to the inverse of the time spent observing the direction 
r (the small-scale details of the pointing of the antenna are irrelevant, as the beam 
smearing will ensure that the function V is smooth on angular scales below the 
beam width). Secondly, let us ignore the nuisance terms®’^ from monopole, dipole, 
etc., so that we can omit A from equation (1) of Einally, we let denote the 
power spectrum that results after the experimental beam smearing has been taken 
into account (this is simply the real power spectrum multiplied by some coefficients 
that approach zero for very large £, corresponding to angular scales far below the 
beam width). 

Let us focus on some fixed value 1 = 1* for now, say 1* = 17. The analysis is 
then repeated from scratch for all other Avalues of interest. Since the power Ci* has 
units of /uK^, we clearly want to estimate it by some quantity Ci* that is quadratic 
in the data. The simplest estimator of this type is 




'ifj(r)x(r)dQ 


( 1 ) 


where ip is some function on the sphere, and the most general estimate is readily 
shown to be just a weighted average of such estimates^. A straightforward calcula¬ 
tion shows that 


OO I p 

{Ci*) = Y. E + / |V’(?)|V(?)dQ, (2) 

.^=0 m=—l 

where ipim denotes the spherical Eourier transform of ip, i.e., the coefficients in an 
expansion of ip in spherical harmonics. In other words, we see that the expectation 
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value of our estimator is the sum of two terms of quite diiferent character. The 
first, the contribution from cosmology, is the power spectrum convolved with a 
window function = J2m l'0/mP- The second, the contribution from noise, is 
just an average value of V, the weights being This is very similar to the 

result when estimating the power spectrum from a galaxy survey®. Just as in that 
paper, we will find it very convenient to use the standard Dirac quantum mechanics 
notation with kets, bras and linear operators. This allows us to write = (^'m.1'0). 
A window function should always integrate to unity, so the correct normalization 
for Ip is just (V’lV’) = 1- Defining the operator 

OO I 

L = ( 3 ) 

.^=0 m=—l 

equation (2) becomes simply 

{Ci*) = {MCL + V{vm. (4) 

Note that T is a scalar operator satisfying L\lm) = l\lm), and is related to the 
(vector) angular momentum operator L = —ir X V through = L{L + 1). 

Now what is the best choice of 'ipl Equation (1) tells us that is the square 
modulus of a random variable whose real and imaginary parts are both Gaussian. 
Thus if Ip is real, the standard deviation of Ci* (the vertical error bar) is simply 
\/2 times its expectation value. (If the real and imaginary parts contribute equally, 
this decreases by a factor of \/2.) Basically, we minimize the vertical error bars by 
minimizing {Ci*). But assuming that our window function is narrow enough that we 
are measuring mostly what we want to measure, Ci, the first term in equation (4) 
satisfies {'iP\Cl\'>P) ~ Ci, independent of 'ip, so we minimize the vertical error bars 
by simply minimizing the second term, {'ip\V(r)\'ip). As a measure of the horizontal 
error bars, we will use A£, the r.m.s. deviation of the window function from £*. 
With our quantum notation, we have simply A£^ = {'ip\{L — £*)^\'ip). As we will 
see, it is impossible to minimize both error bars at the same time, since there is 
a trade-off between them. It would be like asking for the best and cheapest car. 
Instead, the best we can do is minimize some linear combination E = (H), where 
we have defined 

H = (L-£*)^+l3Vir), (5) 

and the parameter /3 specifies how concerned we are about the vertical error bar 
relative to the horizontal one. Continuing our quantum analogy, we see that we 
want to find the 'ip that minimizes the total “energy”, where the “kinetic energy” 
(T — £*)'^ corresponds to the horizontal error bar and the “potential energy” l3V{r) 
corresponds to the vertical error bar. If we for sake of illustration set £* = 1/2, 
we simply want to minimize + (3V{r)\'ip), given the constraint {'ip\'ip) = 1. 

Introducing a Lagrange multiplier E, we arrive at the Schrodinger equation 

[L^+(3V{r)]\i’) = E\i’)- ( 6 ) 

In other words, we want to find the ground state wavefunction for a particle confined 
to a sphere with some potential. From our knowledge of quantum mechanics, we 
can immediately draw a number of conclusions about the solution, all which turn 
out to agree well with the exact numerical results^. 
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• I'i/’p will be small in regions where the noise variance is large, so regions that 
received little observation time will receive low weights in the analysis. 

• Except for the case of complete sky coverage, we will have > 0. 

• If incomplete sky coverage confines •0 to a region of the sky whose angular di¬ 
ameter in the narrowest direction is of order A0, then the uncertainty principle 
tells us that the minimum /S.I must be at least of order 1/A0. 

• This limit on the spectral resolution is independent of 1. 

• For a sky map of COBE type, where A9 is of order a radian given a 20° 
galactic cut, the uncertainty principle thus gives Ai > 1. This agrees well 
with the horizontal error bars actually attained in Figure 1. 

• If the sky-coverage is incomplete, V is infinite outside of the region covered, 
and we recover the quantum-mechanical particle-in-a-box problem. From this 
we know that ip will always go to zero smoothly as it approaches the survey 
boundary. 

• This smoothness of ip is really the gist of the method, as it radically reduces 
“ringing” in Fourier space, “kinetic energy”, without increasing the “potential 
energy” {ip\V{Y)\ip) much at all. 

• If the the survey volume consists of several disconnected parts, then Al is 
limited by the A6 of the largest part. For the galaxy-cut COBE case, for 
instance, using only the northern half of the sky gives the same Al as using 
both the northern and southern skies combined. (However, including both of 
course helps reduce the vertical error bars.) 

What typically happens when using the full optimization machinery^ is that several 
solutions Ip are found to be almost equally good. For the COBE case, for instance, 
the ground state solution for a given i* comes out to be almost degenerate, with 
(21* + 1) orthogonal eigenfunctions all giving almost the same E. This is because 
the azimuthal symmetry of the galactic cut preserves the orthogonality of spherical 
harmonics corresponding to different m-values. When this happens, we obviously 
want to average the (21* + 1) different estimates of C^*, as this reduces the vertical 
error bars by a factor of \/2l* 1 without hardly widening the horizontal ones at 

all. This is basically how Figure 1 was made. 

The author wishes to thank Ted Bunn, George Efstathiou, Carlos Frenk and 
Joseph Silk for useful comments on the manuscript. 
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